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Abstract. We investigate the electronic transport coefficients in unconventional su- 
perconductors at low temperatures, where charge and heat transport are dominated by 
electron scattering from random lattice defects. We discuss the features of the pairing 
symmetry, Fermi surface, and excitation spectrum which are reflected in the low tem- 
perature heat transport. For temperatures fcsT Ss 7 -C -4o, where 7 is the bandwidth 
of impurity induced Andreev states, certain eigenvalues become universal, i.e., inde- 
pendent of the impurity concentration and phase shift. Deep in the superconducting 
phase (fcsT 7) the Wiedemann-Franz law, with Sommerfeld's value of the Lorenz 
number, is recovered. We compare our results for theoretical models of unconventional 
superconductivity in high-T c and heavy fermion superconductors with experiment. 
Our findings show that impurities are a sensitive probe of the low-energy excitation 
spectrum, and that the zero-temperature limit of the transport coefficients provides an 
important test of the order parameter symmetry. 1 

1 Introduction 

Heat and charge transport in metals provides valuable information on the spec- 
trum of charge carriers and phonons, as well as the the scattering of the carri- 
ers by defects and impurities. In a normal metal at low temperatures the con- 
tribution of the phonons to the electrical and thermal conductivities becomes 
negligible compared to the electronic part. When the temperature is further- 
reduced the electronic transport is dominated by scattering of electrons from 
lattice defects. Eventually the electrical conductivity approaches a constant, 
cr(T) oc constant, while the thermal conductivity becomes linear in tempera- 
ture, k(T) oc T. Both transport coefficients are related by the Wiedemann- 
Franz (WF) law k/{<jT) — Ls, with Sommerfeld's result for the Lorenz number 
Ls = ^-(fcs/e) 2 . Inelastic scattering leads to violations of the WF law and to 
significant deviations from L$. 



1 This work was done in collaboration with S.-K. Yip, J. A. Sauls, and D. Rainer. 
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This picture is dramatically altered upon entering the superconducting phase 
with the formation of a pair condensate and the opening of a gap in the exci- 
tation spectrum. The excitation gap in conventional (s-wave) superconductors 
destroys the WF law. This effect has been known since the early days of the 
BCS (Bardeen, Cooper & Schrieffer 1957) theory of superconductivity from the 
work of Mattis and Bardeen (1958) on the electrical conductivity, and from the 
analysis of Geilikman (1958) of the heat transport. 

In this article we investigate the behavior of the heat current for an uncon- 
ventional superconductor, i.e. , for an order parameter with reduced symmetry 
for which gapless excitations exist even at zero temperature. Such supercon- 
ducting states have been argued to exist both in the cuprates and heavy fermion 
superconductors. [For reviews on the heavy fermion system UPt3 see, e.g., Sauls 
(1994a), and Heffner & Norman (1996), and on the cuprates, e.g., Pines (1994a), 
and Scalapino (1995)]. The leading pairing candidate for a tetragonal crystal 
structure (D^) in the cuprates is the B ig state, a spin-singlet state with lines 
of nodes at the Fermi positions pj x = ±Pf y - The most promising candidates in 
the heavy fermion system UPt 3 are the two-dimensional orbital representations 
coupled to a symmetry breaking field. For UPt 3 , which has a hexagonal crys- 
tal structure (Dq^), phase diagram studies (Bruls et al. 1990, Adenwalla et al. 
1990), and transport measurements (Shivaram et al. 1986, Mullcr et al. 1987, 
Broholm et al. 1990, Signorc et al. 1995) lead to either an even-parity, spin 
singlet Ei g state, or an odd-parity, spin triplet E2« pairing state. In both cases 
the order parameter vanishes at the Fermi surface on a line in the basal plane, 
Pf z = 0, and at points at the poles, pf x = pf y = 0. Other models studied in the 
Ginzburg-Landau regime include the one-dimensional orbital, spin triplet model 
with weak spin-orbit coupling by Machida & Ozaki (1991) and Ohmi & Machida 
(1993), the accidental degeneracy AB-model with two unrelated one-dimensional 
orbital representations by Chen & Garg (1993, 1994), and the accidental degen- 
eracy AE-model by Zhitomirsky & Ueda (1996). 

It is known that in a clean superconductor with an order parameter that 
vanishes along a line on the Fermi surface the low-energy density of states is 
linear in the excitation energy, N(e) ~ Nf e/Z\ for e < Z\ , which leads at low 
temperatures to a power law behavior in transport properties as a function of 
temperature for T -C T c (Coffey et al. 1985, Pethick & Pines 1986, Barash et al. 
1996). Gor'kov & Kalugin (1985) and later Choi & Muzikar (1988) showed that 
this spectrum is altered by a random distribution of impurities. A new low energy 
scale, 7, emerges below which the density of states is approximately constant and 
non-zero at zero energy. The energy scale 7 is interpreted as the bandwidth of 
quasiparticle states bound to impurities (Shiba 1968, Rusinov 1969, Buchholtz 
& Zwicknagl 1981, Prcosti et al. 1994, Balatsky et al. 1995). These impurity- 
induced states develop below the superconducting transition. They are formed by 
the constructive interference of particle- and hole-like excitations that undergo 
Andreev scattering from the anisotropy of the order parameter in momentum 
space. For an order parameter with a line of nodes that is required by symmetry, 
the bandwidth and density of Andreev states at zero energy are finite for any 
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finite concentration of impurities, ^ 0. Both 7 and N(0) depend on the 
impurity concentration m and the scattering phase shift 5 . Thus, 7 defines a 
crossover energy scale or crossover temperature T* , below which the transport 
properties are dominated by the Andreev states. For excitation energies above 7 
the transport properties are determined primarily by the scattering of continuum 
excitations. 



We study in detail the low temperature behavior of the thermal conductivity 
tensor for unconventional superconductors with line and point nodes in the order 
parameter. One of the issues we address is the universality of the electrical and 
thermal conductivities at low temperature. We show that the components of the 
thermal conductivity tensor that depend on quasiparticles near the line nodes 
are determined by the same scattering rate as the electrical conductivity and are 
universal in the limit T — > 0. Furthermore, we show that the WF law is recovered 
in the limit T <C T* , independent of the universality of the transport coefficients. 
However, a significant temperature dependence of the Lorenz ratio occurs over 
the temperature range, T* <, T < T c , even for elastic scattering. The universal 
eigenvalues for k and <x result from the cancellation between two factors: (i) 
the density of Andreev states, which is proportional to N(e^> 0) <~ Nf-f/A , 
and (ii) the reduction of phase space for scattering of gapless excitations, which 
is proportional to r(e — > 0) <~ ft/7. This leads to estimates for the electrical 
conductivity, <r(T) ~ N(e) vj r(e) ~ NfVj(U/A ), and the thermal conductiv- 
ity, k(T) ~ N{e)k 2 B TvjT{e) ~ N f vfk 2 B T(h/A ), which arc independent of the 
defect density or scattering phase shift. Perhaps the most surprising result is 
that the ratio of the thermal and electrical conductivity gives the Sommerfcld 
value for the Lorenz ratio, njoT ~ L$ = ^-(fcs/e) 2 . Thus, the differences in 
the coherence factors that determine the conductivity tensors, k and er, do not 
affect the Lorenz ratio L(T) = k/ctT for r«T* and huj <C 7. This discovery 
was made independently by Sun & Maki (1995) for a two-dimensional d x 2_ y 2 su- 
perconductor. For temperatures above the crossover region, T > T*, the Lorenz 
ratio deviates significantly from the Sommerfeld value. The WF law breaks down 
due to the energy dependence of the scattering lifetime in the superconducting 
state, although scattering of electrons from impurities is purely elastic. 

The rest of this article is organized as follows. First we introduce the lin- 
ear response theory for charge and heat currents in unconventional supercon- 
ductors. Then we discuss in general the conditions necessary for observing the 
Wiedemann-Franz law deep in the superconducting state. Next, several order 
parameter models are introduced and the low-temperature behavior of the corre- 
sponding electrical and thermal conductivities are discussed. Finally, we compute 
the thermal conductivities in the superconducting phase and compare our results 
with recent experimental measurements in UPt 3 and Zn doped YBa 2 Cu306.g 
single crystals. 
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2 Electrical and Thermal Conductivities 



We consider a superconductor with anisotropic singlet pairing or unitary triplet 
pairing and discuss the electrical and thermal conductivities in the long wave- 
length limit q^O, and at temperatures T^O. For simplicity we assume isotropic 
impurity scattering, which we treat self-consistently to leading order in l/fc/£o- 
In this case the first order corrections to the current response functions of the 
impurity self-energy, 8<7i mp , and the order parameter, 8 A, vanish for all pairing 
states listed in Table 1. Self-energy corrections corresponding to the excitation 
of collective modes of the order parameter, 8 A, also vanish in the limit q — > 
(Hirschfeld et al. 1989, Yip et al. 1992). The self-energy corrections are also 
called 'vertex corrections' in the Green function formulism of the Kubo response 
function (Rickayzen 1976). If vertex corrections do not contribute, the response 
functions for the electrical and thermal conductivities depend only on the equi- 
librium propagators and self-energies, and the external perturbations. A detailed 
derivation and justification of the linear response in the quasiclassical formula- 
tion of a Fermi liquid has been given earlier by several authors and is repeated 
in brief in the Appendix. Hence we simply present the response functions and 
start with the discussion of our results for various order parameter models. 
For spin singlet states the electrical conductivity is given by 

e 2 N f h f f 

Rfi<ry(w,T) = —^~ jde J dp f v u v fd [f (e_) - / (e+)] 
x Re |m*( P/ ; e, u) (<7*( P/ ; e_) ff *( P/ ; e+) + ( P/ ; e_)/ fl ( P/ ; e+) + tt 2 ) 

-M«( P/ ;e,a;)(^( P/ ;e_)^ (1) 

where e± = e ± Ku>/2, and /(e) is the Fermi-Dirac distribution function, and 
Jdpf ... is a normalized Fermi surface integral. For triplet pairing the spin scalar 
product is replaced by its corresponding spin vector product, /o/o — * f • fo- The 
retarded (R) and anomalous (a) auxiliary functions M x are defined as (for more 
details see the Appendix), 



M x (p f ;e,u)= — - + ' ^ 2 forXeRa}. (2) 



C+(P/;e,") 
^CX{p f -e,uY + D x (e^y 



The result for <Tij(cu,T) was obtained earlier for electron-phonon and impu- 
rity scattering in conventional superconductors in the strong coupling limit by 
Lee et al. (1989) and for the in-plane conductivity of layered weak-coupling 
superconductors by Graf et al. (1995). The conductivity formula reduces for 
a dirty s-wave superconductor to the well-known result of Mattis and Bardeen 
(1958). 
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In contrast to er^- , the thermal conductivity for a spin singlet state is deter- 
mined solely by the anomalous part of the response function, 

x M-Gv ; e, 0) ( ffo A ( P/ ; e)^( P/ ; e) - e)/ * ( P/ ; e) + . (3) 

The retarded and advanced contributions dropped out after applying the normal- 
ization condition. Physically, this means that the deviation of the quasiparticlc 
distribution function due to a thermal gradient determines the heat current; 
changes in the quasiparticle and Cooper pair spectrum do not. Equation (3), 
combined with the equilibrium propagators, impurity self-energy and order pa- 
rameter, is the basic result for the electronic contribution to the thermal con- 
ductivity tensor. It can be shown that Kij in equation (3) reduces to the same 
expression for the thermal conductivity as reported earlier by Schmitt-Rink et 
al. (1986), Hirschfcld et al. (1986, 1988), and Fledderjohann & Hirschfeld (1995), 
except that these authors do not include the D°_ term from the impurity self- 
energy. The terms are only significant at low temperatures, where they can 
help to distinguish between magnetic and nonmagnetic impurities, and vanish 
in both Born and unitarity limits. 



2.1 The Wiedemann- Franz law 

In the limit T — > and lu — » the occupation factors /(e— — ) and 

df(e)/de confine the e-integrals in (1) and (3) to a small e-region of order fc B T 
or hui. Assuming the existence of an energy scale e* » fc^T, on which the 
propagators and self-energies vary, we can set e = in the slowly varying parts 
of the integrands. Using the normalization condition, g^(pf) 2 — /^(P/)/(f(P/) = 
— 7r 2 , in addition to the general symmetry relations for the Green functions, we 
obtain for the conductivities 

= ^J^-M)J dpivl , lv/J ^g- y (4) 

and 

It is useful to write our final results in terms of a mean Fermi velocity and an 
effective transport scattering time r^ , which incorporates all of the coherence 
effects of superconductivity at T — > 0. 

The energy integrals are standard, so the conductivities for a system with D 
dimensions reduce to 

2 

Reaijiu) -> 0,T-> 0) = e —NfVfjTijVfj , (6) 
Kij (T - 0) = ^NfVf^Vfj , (7) 
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where Vj { = Jrfp/V/ ; j(p/) 2 , and the effective transport time is defined by the 
tensor 

Here e fl (e) is the impurity renormalized quasiparticlc excitation spectrum. 

For an isotropic normal metal one has e R (0) — ih/2r, where t is the quasi- 
particlc lifetime due to impurity scattering in the normal state. The transport 
lifetime in the normal state reduces to r for isotropic impurity scattering, i.e., 
Tij — rSij. Note that expression (8) is applicable to the normal state because 
the key assumption in deriving (6) and (7) was that T is small compared with 
e*, where e* is the energy on which the propagators and self-energies vary. Thus, 
for the normal state e* ~ Ef, while for the superconducting state e* ~ 7, where 
7 is the impurity bandwidth of Andreev states, defined by the transcendental 
equation 



nN f cot 2 S + ( 7 / v /|^(P/)l 2 +7 2 ; 2 



with the Fermi surface average (...), the scattering phase shift So, and the im- 
purity concentration n». 

In some respect the band of impurity states forms a new low-temperature 
metallic state deep in the superconducting phase. This analogy becomes more 
transparent when we calculate the temperature corrections to the transport coef- 
ficients, using a Sommcrfcld expansion. However, the 'metallic' band of impurity 
states has other properties that differ significantly from those of conventional 
metals. The special features of the impurity induced metallic band reflect the 
reduced dimensionality of phase space for scattering and the energy dependence 
of the particle-hole coherence factors. Both features combined lead to (i) the 
universality of the transport coefficients at T — > for excitation gaps with line 
nodes or quadratic point nodes, and (ii) the temperature dependence of the 
Lorenz ratio for elastic scattering at T* < T < T c . 

We emphasize that (6) and (7) hold for gapless superconductors in which 
the leading contribution to the transport current is that from quasiparticle ex- 
citations with energies T & T* <C T c . For superconductors with a gap at the 
Fermi surface the number of quasiparticle excitations is of an activated type at 
low temperatures, ~ exp (— Ao/T), and thus the transport coefficients cannot be 
described by equations (6) and (7), with the consequence that the Wiedemann- 
Franz law is strongly violated below T c . 



2.2 Order parameter models 



In our analysis we take a very simple approach. Instead of examining the effects 
of a multi-sheet Fermi surface on the heat and charge current, we model the 
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excitation spectrum by an excitation gap that opens at line and point nodes 
on the Fermi surface, and by the Fermi surface properties in the vicinity of the 
nodes (i.e. the Fermi velocities, v/, and the density of states, Nf, near the nodes). 
The low-temperature behavior of the transport coefficients probes these 'lower- 
dimensional' regions of the Fermi surface, e <C A a <C Ej, where the excitation 
gap vanishes, and is less sensitive to the overall geometry of the Fermi surface. 
It is therefore admissible to parametrize the nodal regions of the gap with a 
minimal set of nodal parameters. The advantage of this approach is that we can 
quantitatively determine the phase space contributing to the low-temperature 
transport coefficients and then examine in more detail the effects of impurity 
scattering and order parameter symmetry on the heat current, without having 
to know the overall shape of the Fermi surface or basis functions. 

We consider models for the order parameter belonging to several irreducible 
representations of the D 6 h and D±h symmetry groups, which are chosen because 
they represent the primary candidates for the pairing symmetry of the heavy 
fermion superconductor UPt3 [for more details on the heat conduction at low 
temperatures for the various order parameter models see Graf et al. (1999)], 
and the high T c superconductor YBa 2 Cu 3 06.9. 2 The singlet and triplet pairing 
states, we consider, have the form 

(singlet) A(p f )=A (T)y r (p f ), (10) 
(triplet) A(p f )=A (T)y r (p f )d, (11) 

where d is the quantization axis along which the pairs have zero spin projection. 3 
The basis functions yp listed in Table 1 have nodes in the order parameter 
and excitation gap that are dictated by symmetry. We model the different nodal 
regions of the gap by a set of parameters, that determine the variation of 
the gap near the nodes (Graf et al. 1996a). Consider the E-rep models shown in 
Fig. 1. In the vicinity of the equatorial line node (pf z = 0), \A(9)\ ~ /iii nc Z\ |6>|, 
while near the poles (pf x =pf y = 0), \A(9)\ ~ Upoint A \O\ n , where the internal 
phase winding number is n = 1 for Ei g and n — 2 for E 2u . 4 

2.3 Universal limits and low temperature corrections 

For simplicity we restrict our discussion of the universal limits and low-T correc- 
tions to the E-rep models proposed for UPt3, and list the asymptotic values for 

2 For detailed discussions of the possible pairing models in UPta and high T c super- 
conductors see the reviews on UPt3 by Sauls (1994a) and Heffher & Norman (1996), 
and on the cuprates by Pines (1994a) and Scalapino (1995). 

3 In this analysis spin has no effect on the transport coefficients beyond the connection 
between the direction of the spin quantization axis and the nodal structure of the 
odd-parity, triplet state basis functions. 

4 The parameters Mi; ne and /i po int define the slope or curvature of the gap near a 
nodal line or point in a spherical coordinate system that is obtained by mapping an 
ellipsoidal Fermi surface onto a sphere. 
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Table 1. Symmetry groups, irreducible representations, and the standard basis func- 
tions of the low temperature phase of several pairing models. 
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the other pairing models in Table 2. For a uniaxial, unconventional supercon- 
ductor with an ellipsoidal Fermi surface a Sommerfeld expansion of the charge 
and heat conductivity tensors gives 

<Ti(T) ~ a CTj j + fi a> i T 2 , i — n h r fl 21 

In the strong scattering limit the parameters a fc and (3 b are related to the mi- 
croscopic model parameters by 



a K , b ~ — 7s rzv , — ~ — — =- , (13) 



3e 2 /3 CT;6 5 /? Ki6 
a <^.6 - ~TTT a «.6 ' ■ I 14 ) 



Here 7s = |7r 2 fc 2 J 7V/ is the normal-state Sommerfeld coefficient, and the effective 
transport scattering time is defined by r a = t\> = ta, where 

t a = 37i/(4/wV0)) . (15) 

For heat flow along the c-axis the coefficients depend sensitively on the inter- 
nal phase winding number, n, of the order parameter around the point nodes, 

a ^ c ~ 3 4^ poin tZ\o(0) t 1, (E 2u ) ' [lb > 

3e 2 

- -5T2-a K ,c , (17) 

7T K B 
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0.0 0.5 1.0 



(71/2) 

Fig. 1. The excitation gap (exaggerated) at the Fermi surface (sphere) for Ei g (a) and 
E2 U (b) order parameters showing the line node in the basal plane and the linear and 
quadratic point nodes at the poles. Panel (c) shows the normalized excitation gap at 
T — for both states as a function of the polar angle O. The parameters correspond 
to those that fit the low temperature thermal conductivity data shown in Fig. 4. 



^1 ~ 7n2k B x / 5 > (Eig) ^£ ~ 1^6 / 3 , (E lg ) 

a K , c - 120 7 2 \2,(E 2u ) ' a a , c ~7a K , b \ 5 , (E 2u ) ' 1 ° J 

One important result is that at low temperatures the basal plane transport mea- 
surements cannot distinguish between the two pairing models, whereas c-axis 
transport coefficients are sensitive to the symmetry of the pairing state. Small 
variations in the concentration of defects may be used to probe the symmetry of 
the pairing state. This is evident in the asymptotic slope of the thermal conduc- 
tivity a K , c = liniT^o K c/T, which is universal for the E 2m state, but depends on 
the impurity concentration (rij) for the Ei g state. Furthermore, in the unitarity 
scattering limit 7 2 oc n^; thus, the coefficient of the T 3 term scales with the 
impurity concentration as (3 C <~ 1/m for the E 2u state and (3 C <~ for the 
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Ei g state. 

The leading order finite temperature corrections to the Wicdcmann-Franz 
ratio become 

(19) 

which increases with temperature for T T* <~ j/ks- This behavior arises 
from two sources: (i) the residual density of states at e = 0, N(0) ~ iVj (7/A0), 
and (ii) the difference in the coherence factors between thermal and electrical 
conduction. Note that for weak scattering, or very clean materials, the very 
low-temperature regime T < T* may be difficult to achieve in practice. 

The asymptotic limits for the order parameter models listed in Table 1 are 
summarized in Table 2. Note the universal result for the in-plane thermal con- 
ductivity of the d x 2_ y 2 model for the high T c superconductors. We omit the 
asymptotic value of the out-of-planc conductivity in the cuprates because it 
strongly depends on the model for the interlayer transport (i.e. coherent vs. 
incoherent propagation of quasiparticles). 



L(T) 



«i(T) 



T Re <7fc (T, w — > 0) 



L s 1 + 



a K a„ 



Table 2. Asymptotic values of the thermal conductivity tensor n/T at T — > 0. 

Pairing State Rep (Group) n b (T) (JpsTvf^ « e (T) {IflsT vj^ 

d x 2_ y 2 Bi g (D 4h ) 2/(7iy*A) 

hybrid-I Ei g (D 6h ) 1/(2 M A) 7/(M?A?) 

hybrid-II E 2u (D 6h ) 1/(2/* A)) l/(2/x 2 A) 

spin triplet B iu (D 6h ) 3/(2/*Ao) l — (//3A0/7)* 

jl*3A) 

AB-model A 2u ffi jBi u (D 6h ) 3/(2/* A) - ^^Zllf f or ^ _ i 



AE-modcl Ai s iEi K (D 6h ) 



Ms A) 

7 3 ~7 



(/*aA?)(MS^) 3 (/^AfX^A?) 



The low temperature anisotropy ratio k c (T)/ Kb(T), suggested as a symmetry 
probe by Fledderjohann and Hirschfeld (1995), distinguishes between the differ- 
ent E-rep pairing models and is very sensitive to the impurity concentration, see 
Table 2. The experimentally observed anisotropy ratio in UPt 3 (Lussier et al. 
1996) approaches the value 



(20) 
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Assuming scattering in the unitarity limit and a reasonable value for the scatter- 
ing rate r(T c ) = H/2t(T c ) w 0.06 . . . 0.12 fcflT c , we can account for the data with 
both pairing states Ei g and E2„ by adjusting the nodal parameters at the poles. 
We also can account for the data within the AE-model of Zhitomirsky & Ueda 
(1996) by fitting the nodal parameters of the crossing line nodes (cross-nodes). 
States with cubic point nodes, e.g. the spin triplet model B lu by Machida and 
coworkers (1991, 1993) and the accidental degeneracy AB-model A2 U © «B lu by 
Chen & Garg (1993, 1994), have a much larger phase space contribution from 
the quasiparticles scattering from impurities near the nodes. Further experiments 
on UPt3 with controlled studies for various impurity concentrations should dis- 
tinguish between the E-rep and AE models as shown for the E-rep models in 
Fig. 2. 
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Fig. 2. Impurity dependence of the thermal conductivity ratio for an Ei a and E211 order 
parameter. The symbols are the low-temperature fits of the impurity concentration for 
the UPt3 single crystal measured by Lussier et al. (1996). 



3 Numerical results 



We explore the temperature dependence of the thermal conductivity for un- 
conventional superconductors. Similar results were obtained earlier by several 
authors (see e.g. Schmitt-Rink et al. 1986, Hirschfcld et al. 1986, 1988, 1994, 
Monien et al. 1987, 1988, Arfi et al. 1988, 1989, Won et al. 1993, and Sun et 
al. 1995). What sets this work apart is the focus on the impurity induced band 
of Andreev states, and their influence on the low-temperature behavior of the 
heat current. Here we extend our analytic analysis with numerical calculations 
of the low-T behavior beyond the universal region. The numerical results are 
obtained by computing the equilibrium propagator, scattering self-energy, and 
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order parameter self-consistently. These results are then used to compute the 
transport coefficients. 

3.1 Thermal conductivity in YBa 2 Cu30 6 . 9 

In Fig. 3 we plot for a d x 2_ y 2-wave pairing state, A(jpf) = A cos 2<f>, the ther- 
mal conductivity for several normalized scattering cross sections a = sin 2 Sq, and 
normalized scattering rates a — h/ (2nkBT c0 T). 5 A consequence of the universal 
limit is that the ratio 



scales linearly with the scattering rate, -T, and is independent of the scattering 
strength a. For T £ T* Arfi et al. (1989) have shown that k(T)/T oc (1 - 
<j)k(T c )/T c strongly depends on the scattering phase shift. This explains (i) 
the sudden drop of [k(T)T c ]/[k(T c )T] in Fig. 3(a) at ultra-low temperatures 
for weak scattering, where the universal limit is achieved only for T < T* ~ 
(2/iii ne Z\o/fcs) exp (— ir[iu ne Ao/4r), and (ii) the scaling of the zero temperature 
intercept with a in Fig. 3(b). Below T* the ratio k/T approaches the universal 
limit. Weak scattering leads to an approximately linear temperature dependence 
over a large portion of the temperature range (Pethick et al. 1986, Barash et al. 
1996). However, k/T changes drastically in rather clean superconductors below 
the exponentially small crossover temperature, where it approaches its linear 
low-temperature asymptote. 

Our calculations of the effects of the impurities on the low-T thermal conduc- 
tivity assume a cylindrical Fermi surface and purely elastic scattering processes. 
These simplifications don't affect the low-temperature thermal conductivity since 
in this limit it is governed by the nodal points. However, we do not expect an 
accurate description of the experimental data at higher temperatures, T ~ T c , 
where inelastic scattering is significant and excitations over the whole Fermi 
surface contribute to the transport properties. 

Recently Taillefer et al. (1997) reported measurements of the heat conduc- 
tivity in Zn doped YBa2(Zn :E Cui_ ;r )306.9. The in-plane thermal conductivity 
was measured at temperatures as low as T c /1000, and the predicted universal 
behavior of k/T was observed for the first time. They found that the asymp- 
totic value of k/T is almost independent of the Zn concentration (0% — 3%), 
with limr^o k/T w 18 mW/K 2 m. They also found that phonons scattering off 
the crystal surfaces cannot be neglected in YBa2Cu30g.9 even at T <~ T c /1000. 
The reported asymptotic value is somewhat larger than our previous estimate 
of k/T rs 14mW/K 2 m (Graf et al. 1996b) which is based on the expression 



5 Here we use simultaneously the notation of the dimensionless parameter a and T — 
h/2r for the scattering rate. 



k{T)T c 2Ti AT 



(21) 



T^0K(T c )T 7TT Winc Z\o(0) 7T/Zl ine A)(0) 




(22) 
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Fig. 3. Thermal conductivity of a 2D d x 2_ y 2-w&ve superconductor. In panel (a): k/T 
for different scattering cross sections, a, at a fixed scattering rate a — 0.01. In panel 
(b): k/T for different elastic scattering rates a in the unitarity limit (a = 1.0). The 
isotropic s-wave BCS result is shown for comparison. 



with a Drude plasma frequency lo p w 1.5 eV, and ta — h/nAo(0), but is well 
within the uncertainties of the known material parameters. For example, if one 
uses specific heat measurements instead of the Drude plasma frequency, with 
typical values for the normal-state specific heat C e \/T = js ~ 270±70 J K~ 2 m~ 3 
(Junod et al. 1990), and the effective Fermi velocity near the nodal regions, 
v f w 110 ± 20km/s, then we obtain n/T w 20mW/K 2 m. 



3.2 Thermal conductivity in UPt3 

The inelastic scattering rate in UPt 3 is comparable to elastic scattering at T c 
and dominated by electron-electron scattering in the normal state for 0.5 K < 
T < 2K (Lussier et al. 1995). Elastic scattering is roughly isotropic (s-wave 
scattering). However, recent resistivity measurements by Kycia et al. (1998) 
show deviations of order 20% from isotropic scattering. In order to describe the 
temperature dependence of the thermal conductivity we model elastic scattering, 
To, and inelastic scattering, by a phenomenological scattering rate T(T) = 
H/2t(T) = T (l + T 2 /T C 2 ), where we neglect the effects of anisotropic scattering. 
The inelastic scattering rate decreases below T c even faster than T 2 because of 
the onset of superconductivity and the opening of a gap over most of the Fermi 
surface. We will neglect this reduction in T(T), since we are mostly interested 
in the low-temperature regime, e.g. for temperatures T ^ T c /10 the inelastic 
scattering rate is less than 1% and can safely be neglected. However, we do 
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not expect an accurate description of the experimental data in the temperature 
regime, T c /2 ^ T < T c , where inelastic scattering is not negligible. 

Our analysis of the universal regime implies rather stringent restrictions on 
an equatorial line node 

l im o?7^i"l / ( I C lnr *BT**0.2vWlo(0)ro. (23) 

T^O K b {l c )l 2^i ine Zio(Oj 

The experimentally imposed constraints are consistent with a nodal parameter 
Mine = 2, an intercept of [n b {T) T c ]/[K b (T c ) T] <, 0.02 at T — > 0, and a crossover 
temperature T* <, T c /12 (Lussier et al. 1996, Huxley et al. 1996, Suderow et 
al. 1997). The elastic scattering rate r is adjusted for a best fit to the low-T 
region for the basal plane thermal conductivity. For an E 2m pairing state a good 
fit is obtained for _T = 0.0l7rT cO ~ 0.03 T c , for which the theoretically computed 
intercept is 

t^o n b {T c )/T c y ' 

This scattering rate fits the low-temperature part of the thermal conductivity in 
the superconducting state, but is smaller than that estimated from the normal- 
state transport data and de Haas- van Alphen data, i~o ~ 0.1 T c — 0.2 T c , (Taillcfcr 
et al. 1987). Larger discrepancies were reported earlier by Lussier et al. (1995, 
1996), and Flcdderjohann et al. (1995). Finally, the nodal parameters defining 
the excitation gap near the point nodes are adjusted to fit the heat current along 
the c-axis, n c (T)/T at T <C T c . These values are \i\ = 2/3 for the Ei s model 
and fj,2 = 4 for the Fi2u model. 

A similar analysis for two of the other proposed pairing models for UPt 3 is 
shown in Fig. 5. It is in principal possible to fit the basal plane thermal conduc- 
tivity, however, a fit for the heat current along the c-axis generally fails because 
of the large phase space contribution from the cubic point nodes and cross-nodes 
(crossing line nodes at the poles) for the B\ u and A 2u © iB\ u order parameters. 
This is the case for all order parameter models other than the E-rep models and 
the AE-model. Thus, the only models that are capable of explaining the existing 
data for the anisotropy and temperature dependence of the heat conductivity in 
UPt3 are the 2D E-rep models and the AE-model. 

The existing low-temperature data for the thermal conductivities of UPt3 
(Lussier et al. 1996) are equally well described by either an Ei ff or E 2m or 
AE model. However one possibility to distinguish between these models is by 
studying the anisotropy ratio k c / n b in the low-temperature region for different 
scattering rates. The corresponding zero temperature limit for the E 2u (hybrid- 
II) state is universal, i.e. independent of impurity scattering for 7 <C Aq. It 
depends only on the nodal (material) parameters through the ratio /Uiine/Mpoint- 
In contrast the anisotropy ratio for the Ei s state and AE pairing state depends 
sensitively on the scattering rate and scattering strength, see Fig. 6. 

We summarize the quality of our fits of the thermal conductivity and the 
observable differences between the various models for the pairing state in UPt3 




Fig. 4. Heat flow along the b and c axes in the unitarity limit. The experimental 
data (symbols) were taken from Lussier et al. (1996) and normalized to the values 
at T c ~ 0.5 K. The insets show the normalized anisotropy ratios. Panel (a): Ei 9 -state 
for two scattering rates r(T c ). The slope parameters at the linear point nodes and 
at the line node were fixed to /ii = 2/3 and /ii = 2.0. Panel (b): E2 U -state with a 
scattering rate i~b = O.Ol7rfeBT c o and different nodal parameters at the point nodes, 
[ii = 4.0, 4.5, 5.0. The slope of the line node was fixed to fj,\ — 2.0. 



by plotting the anisotropy ratio K(T) = [n c {T) / n h {T)]/[n c {T c ) / n b {T c )]. Figure 7 
shows that the ratio 1Z(T) for the data of Lussier et al. (1996) is in agreement 
with the pairing states for the Ei a , E 2tl , and AE models. The spin triplet model 
Bi„ and the AB-model cannot account for the anisotropy at low temperatures, 
1Z < 1. Independent measurements of the heat currents by Suderow et al. (1997) 
differ to some extend from those by Lussier et al. (1996), particularly in the 
ratio 1Z. This discrepancy is mainly due to differences in the heat current along 
the c-axis. This may be caused by crystal imperfections along the crystal c-axis, 
and needs further investigation. 



4 Conclusion 



At very low temperatures heat transport is dominated by elastic scattering of 
quasiparticlcs in lower dimensional regions of the Fermi surface in the vicin- 
ity of the nodes of the order parameter. As a result the low energy spectrum 
can be determined accurately without complete knowledge of the Fermi sur- 
face. The residual density of states at zero energy and zero temperature leads 
to a Wiedemann-Franz law for electronic transport deep in the superconduct- 
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Fig. 5. The thermal conductivity for heat flow along the b and c axes in the resonant 
impurity scattering limit. The material parameters are the scattering rate, Jo, and 
the nodal parameters, /ii, at the point (p) and line (I) nodes. The theoritcal results 
are compared with the experimental data (symbols) for UPt3. Panel (a): fit for an 
odd-parity Bi u -state. Panel (b): fit for an accidentally degenerate A 2u © iSi u -state. 




Fig. 6. Thermal conductivity ratio for a set of impurity scattering rates 7~b and a fixed 
r in = 0.0l7rfc s T c0 (T/T c ) 2 for an Ei 9 (a) and E 2u (b) pairing state. 
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Fig. 7. The anisotropy ratio of the thermal conductivities for various 'best fit' pairing 
states. The line for the polar state Ai u is shown for comparison. The experimental data 
are renormalized relative to the upper transition (T c + ). 

ing phase. The results for the low-temperature transport coefficients can shed 
new light on the strength of the scatterering centers because of the sensitivity 
of the coherence factors to the scattering phase shift. Furthermore, our theo- 
retical results suggest that measurements of the heat transport along different 
crystal axes can differentiate between the various pairing models that have been 
proposed for heavy fermion and other unconventional superconductors. 

Excellent agreement with experiments on UPt3 is obtained for pairing states 
with Ei g or E 2tl symmetry at low temperatures, as well as for the AE-model, 
whereas no agreement (along the c-axis) could be obtained for the paring states 
with Biu and A 2u © iBi u symmetries. The differences in the excitation spectrum 
for Ei g and E 2u can be observed at very low temperatures, T T* ~ T c /12, 
when the impurity-induced Andreev band determines the heat conductivity. 
Measurements of the thermal conductivity ratio, K c /nb, at very low tempera- 
tures (T <C T c ) for various impurity concentrations should further differentiate 
between these models; the E 2u model is predicted to have a universal ratio, while 
the ratio for the Ei g state depends strongly on the concentration and scatter- 
ing strength of impurities. At high temperatures, T <~ T c , our calculations show 
deviations from the experimental data, which reflect our neglect of (i) the effect 
of the onset of superconductivity on the inelastic scattering rate, and (ii) the 
splitting of the superconducting transition temperature. Finally, we note that 
our prediction of a universal n/T at low temperatures is in good agreement with 
recent experiments on Zn doped YBa 2 Cu 3 06.9 (Taillefer et al. 1997). 
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A Linear Response in the Quasiclassical Formulism 
A.l Quasiclassical transport coefficients 

The electrical current density is determined by the Keldysh propagator, the 
Fermi velocity, v/ , and the density of states per spin, Nf, at the Fermi level, 

j e (R,t) = 2N f J dpfj A( eV/(Rf )) iTrf3g*( P/ ,R;e,t), (25) 

and the heat current density has the form 

j e (R,t) - 2N f J dp f J A ( £V/ ( P/ )) ±Trg K (p f ,R;e,t) . (26) 

For weak disturbances from equilibrium the current response is linear in the 
applied perturbation. Here, we are interested in the electrical and heat current 
response functions, which define the conductivity tensors by 

5j e = a • E w , and 5j e = -k ■ VT . (27) 

We develop the quasiclassical linear response equations for the transport coeffi- 
cients k and er, following the discussion by Rainer and Sauls (1995), and Graf 
et al. (1995). 

The advanced, retarded, and Keldysh propagators are calculated from the 
quasiclassical transport equations 

[eT 3 -a ext ~a , g J o = v f ■ - Vg , (28) 

and 

(er 3 - a ext - a R ) o g K - g K o (ef 3 - a ext - a A ) 

-a K og A + gKoa K = v f ^Vg K , (29) 

which are supplemented by the normalization conditions (Eilenberger 1968, 
Larkin & Ovchinnikov 1968), 

§ R,A Q ~R,A = _ 7r 2 i ) and g* gK + g K o^ = . (30) 
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To linear order in the perturbation the quasiparticles (with charge e) respond 
to an electric field or temperature gradient through the self-energy terms, 

e h 
crext =--v/-Af 3 , or Vf-VgX, for X e {R, A, K} , (31) 

where A(q, uo) is the vector potential describing the transverse field E w = 
(iw/c)A, and g x are the equilibrium propagators. In order to calculate the 
conductivity we solve the transport equations to linear order in the perturbing 
fields. For the heat conductivity the perturbation is the temperature gradient 
VT. Generally, deviations arise from local equilibrium cither from changes in 
the equilibrium propagators Vg RA or from the distribution function V<?o with 

#o(R) = 1 - 2/(e; T(R)) = tanh ( 2fcj ^ (R) ) • (32) 

We consider only spatially homogeneous superconducting states which are in 
equilibrium, and are "unitary", i.e., the equilibrium mean-field order parameter 
satisfies Z\(p/) 2 = — |Z\(p/)| 2 l, where |Z\| 2 stands for either the spin scalar 
product A A or the spin vector product A ■ A. 



A. 2 Solutions of the linearized quasiclassical transport equations 

The deviations of the propagators from their local equilibrium values, 8g x 
9 

equations for the retarded and advanced propagators 



A g x , and 5a x = a x — a x with X € {R, A, K}, satisfy the following linearized 



8g R > A , h R ' A 



idg 



R,A 



9 , o- ext + da 



(33) 



and for the anomalous propagator, 

h R o 5g a - 8g a oh A = (id<P ) o g A - g R o (id<P ) + 8& a o g A - g R o 5& a 

[0 O , a ext } , (34) 



[a ext , <P ] o og A 



9$o 



where h R A = eT 3 ~a R ' A , and d = TiVf V. We introduced Eliashberg's anomalous 
propagator, Sg a , and self-energy, Sd- a , defined by 

5g K = 8g R o <2> - <? o Sg A + 8g a , (35) 
5& K = 5a R o <Z> - <P Q o Sa A + 5a a . (36) 

The transport equations for $g R ' A < a are simplified by inserting the local equilib- 
rium propagators 

R,A 

( h 4 _ , W i th C R ' A = - i y/\AW\* - (~e R > A r , (37) 
where e R ' A and A R ' A determine the equilibrium self-energy, 



~R,A 

9o 



e R > A r 3 



~R,A 



R,A 



(e-i R ' A )h + A +D R ^ A i. 



(38) 
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Finally, the equations for 6g R ' A and Sg a are solved by using the normalization 
conditions 

9o A ° Sg*' A + 6g R ' A o g^' A = , and o <5.r + ^° o ^ = , (39) 

and h R,A in order to move the propagators Sg R ' A and Sg a to the right on the 
left hand side of (33) and (34). The solutions arc 

5g RA = (cf A g«- A + D R - A ) ~\ (-idg*' A + [a ext + 5a R ' A , g*' A ] J , (40) 
5g a = (C a + g R + D») _1 o ^(iftP ) o g A - g R o (iftp ) 

+(5<r a og A -g R o 5a a - [a ext , <P ] o o g A - g R o [<P , a ea;t ] ^ , (41) 

with the auxiliary functions 

C^( P/ ; £, W ) = C^( P/ ; c+ ) + C RA (p f ; e_) , (42) 

D^( P/ ; e, W ) = D*^( P/ ; e + ) - D RA (p f ; e_) , (43) 

C;( P/ ; e.w) = C*( P/ ; e+) + C A (p f ; e_) , (44) 

I>» ( P/ ; e, w) - D R (p f ; e+) - D A ( P/ ; e_) , (45) 

and with e± = e±Uw/2. Using the normalization conditions (30) we can imme- 
diately invert the matrix 

K^r=-§i^f (46) 

After combining the previous results (40) and (41), according to (35), we 
obtain the general expression for the linearized Keldysh propagator Sg K : 

Sg K = {C R g R + D R ) _1 o + [<7 ext + 6a R , $*] ^ o <2> 

-<2>o o (C^.9 A + D A ) ^ o + [<x ext + ^ , 5 A ] o ) 

+ (C^ + ) o ^(iftP ) o g A - g R o (iftP ) + 5a a o ^ - g R o 6<x a 
- [*ext ,Mo°9 A -90° [*0 , CText] J • (47) 



A closer look at (47) shows that changes in the equilibrium propagators 
dg^' A do not contribute to a static heat current, because the normalization 

condition (30) gives 2 Tv g R ' A dc/Q ' ,A — d Tr (%' A ) 2 = 0. In linear response the 
physical observablcs arc then determined by the deviations from the quasiparticle 
distribution through 8g K . In general the corrections to the self-energies 5a must 
be included and require a self-consistent calculation. 
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